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Particle-wave duality enables the construction of interferometers for matter waves, 
which may complement lasers in precision measurement devices such as gravimeters or 
gyroscopes. This requires the development of atom-optics analogs to beam splitters, 
phase shifters, and recombiners. Integrating these elements into a single device has 
been a long-standing goal. Here we demonstrate a full Mach-Zehnder sequence with 
trapped Bose-Einstein condensates (BECs) confined on an atom chip. Particle inter- 
actions in our BEC matter waves lead to a non-linearity, absent in photon optics. We 
exploit it to generate a non-classical state having reduced number fluctuations inside 
the interferometer. Making use of spatially separated wave packets, a controlled phase 
shift is applied and read out by a non-adiabatic matter-wave recombiner. We demon- 
strate coherence times a factor of three beyond what is expected for coherent states, 
highlighting the potential of entanglement as a resource for metrology. Our results 
pave the way towards integrated quantum-enhanced matter-wave sensors. 



Atom interferometers [JJ are powerful tools for probing 
weak interactions [2] , measuring atomic properties [3] , for 
studying quantum degenerate systems [U [5] , for precision 
measurements [B], and are very sensitive probes for iner- 
tial effects [7J, as beautifully demonstrated for example 
in gravimeters [5] and gyroscopes [§]• 

In most atom interferometers to date, the matter waves 
propagate freely and the two interfering arms correspond 
to wave packets having different momenta |TJ . The beam 
splitters creating the superposition either couple different 
internal states, like in the Raman beam splitter, or leave 
the matter wave in the same atomic state, like in Bragg 
scattering or diffraction from a material grating. 

Bose-Einstein condensates (BECs) are promising can- 
didates for atom interferometry owing to their macro- 
scopic coherence properties [T0~HT3"] . Following the first 
observation of BEC interference in 1997 [TU], various 
building blocks of BEC interferometers have been real- 
ized individually. Interference experiments with BECs 
were performed using Bragg beams in ballistic expan- 
sion [TTJ, with atoms freely propagating in a guide [TU 
115] . and very recently under microgravity [16) . Split- 
ting a single trapped BEC into two separated clouds in 
a double well, interference was observed using ballistic 
expansion after switching off the trapping potential |17l - 
|2"T] . Realizing a full Mach-Zehnder interferometer with 
spatially separated modes of trapped BECs has been a 
long-standing goal 22 22] as the tight confinement pro- 
vided by chip traps allows local probing. 

A fundamental difference between photon and matter- 
wave optics is the presence of atom-atom interactions. In 
dense, confined geometries they lead to dephasing (phase 
diffusion) which ultimately limits the coherence time of 
the interferometers [2TJ |2"5H2"7] . Phase diffusion can be re- 
duced by controlling the interactions with a Feshbach res- 
onance. Coherence times up to 10 s were recently demon- 
strated in Bloch oscillation experiments [35] [5pJ . 



Interactions can also be used to generate non-classical 
correlations between the two paths of the interferometer, 
which can be used either to improve the sensitivity [30 or 
to reduce the effect of phase diffusion |19j . Generation of 
entanglement was recently reported for BECs in multiple- 
well potentials |31] and between two or three internal 
states of spinor BECs 32 36 . In Ref. [32] it was used 
on an internal state BEC interferometer to improve its 
sensitivity below the standard quantum limit <50sql = 
1/y/N, where N is the number of atoms involved in the 
measurement. 



Here we demonstrate a full Mach-Zehnder interferom- 
eter for Bose-Einstein condensates integrated on an atom 
chip. Our interferometric scheme relies on the coherent 
splitting and recombination of a Bose-Einstein conden- 
sate in a tunable magnetic double-well potential, where 
the matter wave is confined at all times. Thanks to a spa- 
tial separation of ~ 2 fjxa between the two wave packets, 
our geometry is sensitive to accelerations and rotations. 
By tilting the double well out of the horizontal plane for 
a variable time t$, we apply an energy difference e and 
thereby imprint a controlled relative phase between the 
interferometer arms. A non-adiabatic recombiner trans- 
lates the relative phase into an atom number difference, 
which is directly read out using a highly sensitive time- 
of-flight fluorescence detector [37] . We show that the use 
of a non-classical state with reduced number fluctuations 
inside the interferometer helps increasing the interroga- 
tion time, which is still limited by interaction-induced 
dephasing. This represents an unambiguous demonstra- 
tion of the interplay between number fluctuations and 
phase diffusion in the case of spatially separated modes 
of Bose-Einstein condensates. 
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RESULTS 

The atom chip used to manipulate the atomic wave 
packets is depicted in Fig. Ek. The interferometer se- 
quence is illustrated in Fig. and Fig. []J displays the 
resulting oscillation of the mean of the normalized pop- 
ulation difference z = (N-^ — Nn)/N t as a function of the 
phase accumulation time t^, where Ni, and Nr are the 
populations of the left and right wells respectively and 
N t = N L + N R . 

In the following we will discuss each step of the se- 
quence in more detail. We first focus on the BEC 
splitter and find that the state generated inside the in- 
terferometer is characterized by a spin-squeezing factor 
f§ = -7.8±0.8dB [3D], implying that it is entangled [38]. 
We then turn to the phase accumulation stage and anal- 
yse the phase diffusion process. We find a threefold ex- 
tension of the coherence time using a number-squeezed 
state inside the interferometer as compared to a coherent 
state. Finally, we introduce the non-adiabatic recom- 
biner and discuss the obtained sensitivity and possible 
improvements. 

Matter-wave source 

The experimental sequence starts with a degenerate 
gas of 87 Rb atoms in the F = 1 hyperfine state in an 
elongated magnetic trap created by our atom chip set- 
up (cf. Supplementary Method 1) [39 . Using radio- 
frequency (rf) evaporative cooling, we prepare samples 
containing N = (N t ) = 1200 atoms at a temperature of 
T ~ 20 nK = h/kg x 420 Hz and a chemical potential 
of fj,/h ~ 500 Hz, &b being Boltzmann's constant and h 
Planck's constant. The brackets denote ensemble aver- 
aging. Given the elongated trap geometry, the cloud is 
in the one-dimensional quasi-condensate regime (cf. Sup- 
plementary Discussion 1) [ID]. 

Splitting 

The entire interferometer sequence relies on rf dress- 
ing [TS1 E5 to dynamically turn the static magnetic 
trap into a double-well potential, whose spacing, bar- 
rier height, and tilt can be controlled by changing the 
amplitude and orientation of the rf field. Splitting is 
performed by increasing the rf amplitude in the dressing 
wires linearly within 5 ms to transform the single trap 
into a symmetric double- well potential [T8J EI]- Along 
the horizontal splitting direction, the potential has two 
minima separated by 2 /Km, the barrier height is 3.7kHz 
and tunnelling is negligible (cf. Supplementary Method 
2). 

To characterize the quantum state obtained after split- 
ting, we measure the distributions of the two observables 



atom number difference n = — Ar and relative phase 
4> = 4>-r. — 4>l- The interferometric sequence can be inter- 
rupted at any time to measure these quantities using de- 
structive time-of- flight imaging (see Methods). The rel- 
ative phase obtained from the Fourier transform of the 
interference fringes is used to monitor the system dur- 
ing the sequence and must not be confused with the one 
measured after recombination (measuring the well's pop- 
ulations). Each experimental run gives access to either 
one or the other observable. We repeat the measurement 
r~> 200 times to obtain good estimates of both the mean 
and the variance of n and <f>. 

For non-interacting particles, the atom number dif- 
ference n is expected to follow a binomial distribution 
with mean (n) — and standard deviation An = y/N 
after splitting. As repulsive interactions energetically 
favour equal populations, reduced fluctuations are ex- 
pected [HI H2] 03J. Figure [2k shows the distribution 
of detected photon number difference s = sl — sr be- 
tween the left and right clouds (our fluorescence de- 
tector registers on average 16 photons per atom, sec 
Methods). Indeed, we observe a suppression of fluctu- 
ations by a factor of two compared to shot noise. After 
correction of the detection noise [371 |H] we deduce a 
number squeezing factor £n = An/^/N = 0.41 ± 0.04 
= -7.8 ± 0.8 dB). This shows that the splitting pro- 
cess yields a number-squeezed state with reduced fluc- 
tuations along one quadrature compared to a classical 
coherent state. 

We now turn to the phase distribution shown in 
Fig. (2Jd. Since the two variables n and <f> can be seen 
as quantum mechanically conjugate {[<f>, n] = 2i) [30] the 
Heisenberg uncertainty relation imposes A(j) > 1/An = 
0.074 rad. We observe a Gaussian-shaped phase distri- 
bution with vanishing average and circular standard de- 
viation g5] A(j> = 0.16 ± 0.01 rad (An A(j> = 2.3), twice 
as large as expected from the Heisenberg lower bound. 
Note that contrary to the number distribution data, the 
phase distribution data were not corrected for imaging 
noise. We conjecture that this broadening is due to finite 
initial temperature [46] and detection noise. 

Altogether, the output state of the beam splitter 
exhibits reduced number fluctuations and high coher- 
ence properties, characterized by the coherence factor 
(cos0) = 0.987 ± 0.001 gTJ. It has been shown gSJ 
that using such a state in an alternative interferomet- 
ric scheme can allow for a sensitivity better than the 
standard quantum limit. The improvement in sensitiv- 
ity is measured by the spin squeezing factor 30 £g = 
£n/(cos0) = 0.41 ± 0.04 and corresponds to a gain of 
7.8 ± 1 dB over the standard quantum limit. Note that 
the numerical value of £s is here identical to that of 
£n thanks to the large coherence factor of almost unity. 
Their statistical uncertainty also coincide because that 
of the coherence factor is negligible compared to that of 
£n- Since spin squeezing, defined as £s < 1 is an en- 
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Figure 1. Scheme of the Mach-Zehnder interferometer and measured output signal, a, Schematic of the atom chip: 
a DC current of 1 A is sent through the trap wire, which, together with a uniform external "bias" field B x — 29.5 G, creates 
an elongated quadrupole trap 60 (jm below the chip. The longitudinal confinement is realized by two wires parallel to the x 
axis (not shown here). An external field -Bo completes the Ioffe-Pritchard configuration. Rf-currents with a relative phase of 
7r are applied on the dressing wires to perform the splitting along x. b, The condensate is coherently split by transforming a 
single trap into a double-well potential; a relative phase between the two arms is imprinted by tilting the double well during 
a time t$; the spacing between the two wells is then abruptly reduced and the potential barrier acts as a beam splitter for 
both wave packets, transforming the relative phase into a population imbalance. After the recombination time Jbs, the atom 
clouds are separated and the particle number in each well is read out by fluorescence imaging, c, The normalized population 
difference between the two wells z = n/N t is measured as a function of t$. It exhibits interference fringes and a damping due 
to phase diffusion. Grey dots: imbalance of individual experimental realizations; black dots: ensemble average (z); red curve: 
theoretical prediction taking into account phase diffusion; dashed black line: expected signal for a classical coherent state for 
which £n — 1- The error bars indicate ±1 standard error of the mean. 



tanglement witness [38] . this also demonstrates that our 
state is non-separable after splitting. The measured fluc- 
tuations and coherence imply that our ensemble contains 
150 entangled particles in the sense of Rcfs. [49, 50 and 
we can exclude entanglement of less than 67 atoms with 
90 % probability. 

Phase shifter 

We now discuss the phase accumulation stage. Af- 
ter splitting, the coherent superposition of left and right 



modes has on average ((f)) — and (n) = 0. We induce 
a deterministic shift of the relative phase by tilting the 
double- well trap out of the horizontal plane while keeping 
the well spacing constant. The tilt is performed adiabat- 
ically with respect to the transverse motion in the trap 
(see Supplementary Note 1) by modifying the relative 
current in two rf wires on our atom chip within 3 ms. It 
is reversed before the recombining sequence. The energy 
difference e [21], which arises from the tilt (see Supple- 
mentary Discussion 2), triggers a linear evolution of the 
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Figure 2. Number and phase distributions after split- 
ting, a, Histogram of the difference between the fluorescence 
signals of the left and right clouds s — sl — sr., in units of 
the standard deviation expected for a coherent state As co h- 
The curves indicate a normal distribution corresponding to 
the measured number squeezing factor £n = 0.41 ±0.04 (solid 
black); the distribution expected in the limit £n = 0, where 
only detection noise is responsible for fluctuations (dotted 
red); and the distribution expected for a coherent state in 
the absence of detection noise (dashed blue) . The inset shows 
a typical fluorescence picture and the regions used to define 
sl and sr. b, Histogram of the measured relative phases cf> 
in units of the circular standard deviation of a coherent state 
A(j> coh . The curves indicate a normal distribution with the 
measured standard deviation A<j> = 5.4 ± 0.5 x A<jf> coh (solid 
black); and the distributions expected for a coherent state 
in the absence (dashed blue) and in the presence (dash-dot 
green) of detection noise. The inset shows a typical matter- 
wave interference pattern from which the phase is extracted. 



relative phase given by 



= ipa + et<j,/h. 



(1) 



In Fig. [3^,, the phase evolution is measured for 7 different 
tilt angles ranging approximately from —9° to +12°. For 
these measurements the phase is obtained directly from 
interference patterns measured after time-of-flight detec- 
tion (cf. inset of Fig.^). Altogether this phase accumu- 
lation stage acts as a tunable phase shifter for the atomic 
interferometer. For the fringes presented in Fig. [TJ; the 
angle was fixed to yield the rate e/h = —349 ± 2 Hz and 
tj, was varied to tune the phase. 

In Fig. [5b-e we observe that the linear evolution of the 



mean phase discussed above is accompanied by a broad- 
ening of the phase distribution, measured by the circu- 
lar standard deviation of the phase A(j> [45] . This phe- 
nomenon illustrates phase diffusion |25j . It corresponds 
to a dephasing arising from the interactions between the 
atoms: since the chemical potential of each condensate 
depends on its atom number, relative number fluctua- 
tions translate into fluctuations of the phase accumula- 
tion rate (cf. Supplementary Note 2 for the influence of 
a residual coupling between the condensates). This leads 
to an increase of the phase variance given by [23 HZ1 IH] 



A^m) = A<t>i+R 2 m 

where the phase diffusion rate R is 



R = 



dN 



N=N/2 



(2) 



(3) 



n(M) being the chemical potential of one BEC with Af 
atoms. The term An = ^y/N corresponds to the stan- 
dard deviation of the atom number difference n. One sees 
from Eqs. (|3| and ^ that number squeezing, defined as 
£n < 1, reduces phase diffusion and allows longer inter- 
rogation times. From a fit (red line in Fig. [3]b) we obtain 
the phase diffusion rate R = 51 ± 4mrad • ms -1 and an 
initial phase spread A0o = 300 ± 90 mrad. The constant 
ti = — 6 ms accounts for the 3 ms used before and after the 
phase accumulation stage to incline and level the double 
well. 

The measured rate is in fair agreement with the ex- 
pected value of R = 46 ± 4 mrad • ms" 1 evaluated from 
Eq. ([3]). Because the splitting process is not adia- 
batic with respect to the axial direction, a slow axial 
quadrupole mode is excited and the condensates are not 
strictly at equilibrium during phase accumulation. Con- 
sequently, the term <9/u/cW in Eq. ^ had to be slightly 
modified to evaluate Eq. (|3| (cf . Supplementary Note 3) . 
This theoretical prediction with no free parameters is rep- 
resented as a grey shaded area in Fig. |3j}. To further 
confirm our understanding of phase diffusion, we directly 
measured <9/i/cW by splitting the initial BEC with an 
angle leading to a controlled population imbalance, lev- 
elling the double well back to horizontal (e = 0), and 
monitoring the evolution of the phase due to a difference 
of chemical potential. Together with the squeezing fac- 
tor and the atom number measured, this gives a third 
value of R = 57 ± 5 mrad • ms -1 also in agreement with 
the two others. For comparison, we evaluate the hypo- 
thetical value R — 112 mrad • ms -1 in the absence of 
number squeezing (£n = 1), indicated as a solid black 
line in Fig. [3J5. The difference with our measurements 
underlines that the interrogation time of our interfer- 
ometer is extended by the fact that the initial state is 
number squeezed. Enhanced coherence times have been 
reported in Ref. [19] and attributed to number squeezing. 
Our measurement of both number and phase fluctuations 
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Figure 3. Evolution of the relative phase and its fluctuations during the phase accumulation stage. All data on 
this figure have been obtained by time-of-flight recombination, a, Linear evolution of the phase for various energy differences 
e induced by tuning the angle between the two wells. The values of e shown in the legend are obtained from linear fits to 
the data (black dashed lines). For all fits, the standard error of the slope is below 3 Hz. The orange points correspond to 
the rate e/h — —349 ± 2 Hz used to record the fringes of Fig. b, Evolution of the circular standard deviation of the 
phase A(j> corresponding to the orange curve of panel a. It exhibits phase diffusion: at short times t<f, < 30 ms the mean 
phase is well defined (blue points); on the contrary, for longer times the phase distribution cannot be distinguished from a 
random distribution (black points, see Methods). The red line is a fit to the blue points with the model of Eq. |2j. Shaded 
area: theoretical prediction without free parameter, taking into account the measured number squeezing (see text). Black line: 
expected behaviour if the initial state were classical (i.e. not number-squeezed). Note that at = phase diffusion has already 
started (see Eq. |2|). c,d,e, Measured phase distributions for three values of indicated by the dashed lines. The error bars of 
a indicate ±1 circular standard deviation divided by \fki, and those of b correspond to ±1 circular standard deviation divided 
by \/2(ki — 1), ki being the number of repetition used for each point i. 



corroborate the direct link between phase diffusion and 
number squeezing. 



Recombiner 

After the condensates have accumulated a relative 
phase, the last step consists in recombining the left and 
right modes in a way that converts the phase into an 
atom number difference. The recombiner plays a role 
similar to that of the output beam slitter of an optical 
Mach-Zehnder interferometer. In a two-mode picture, 
this operation corresponds to a n/2 pulse. In Ref. [2U] 
the phase-dependent heating caused by the merging of 



two condensates was used to infer the relative phase. Us- 
ing Josephson oscillations [51] was also proposed for this 
purpose [H] . Nevertheless one expects the contrast 
of a Josephson recombiner to be fundamentally limited 
by the onset of self-trapping when the interaction energy 
is not negligible compared to the tunnelling energy |52) . 
Here we introduce a different method which makes use 
of a non-adiabatic transformation of the potential. We 
abruptly decrease the rf-dressing intensity to reduce the 
well spacing to 1.5 /zm and the barrier height to about 
1 kHz (see Fig. [4^,-b) . We estimate the energy gap be- 
tween the two lower-lying single-particle eigenstates of 
the double- well potential to be approximately h x 80 Hz. 
The clouds are accelerated towards the barrier and, af- 
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Figure 4. Principle of the BEC recombiner. a,b, Nu- 
merical simulation of the evolution of the transverse in situ 
density when the double-well spacing is abruptly decreased 
(0 < tBS < 2.25 ms) and increased again (tss > 2.25 ms) for 
two different initial relative phases: cj> = (a), and <j) = 7r /2 
(b). Black lines: positions of the two potential minima. c,d, 
This phase sensitive recombination can be understood as the 
propagation and interference of two wave packets launched 
onto a semi-reflective tunnel barrier. The final wave pack- 
ets (filled curves with outgoing arrows) are the result of the 
constructive or destructive interference between the reflected 
(dashed lines) and transmitted (dotted lines) components 
arising from the incoming wave packets. 



ter an adjustable time ^bSj the barrier is raised again to 
separate the atoms for counting. Starting with a state 
having a phase close to 7r/2, we adjusted ies m order to 
maximize the final population imbalance. This yielded 
ies = 2.25 ms. Repeating this operation for different ini- 
tial phases, we observe the sine-like dependence of the 
population imbalance shown in Fig. [TJ;. 

This effect can be understood as the splitting of the 
wave packets on a semi-reflective potential barrier and 
the interferences between them (see Fig. |4j;-d) . Indeed, 
in the absence of interactions and assuming two symmet- 
rically incoming wave packets with relative phase <f>, the 
linearity of the Schrodinger equation and symmetry of 
the potential assures that the population difference be 
given at any time by 



n{t 



BS) 



N t C(t BS ) sin cf>. 



(4) 



The contrast of the read-out — 1 < C(iBs) < 1 depends 
on the time spent in the coupled double well and is a 
function of the overlap of the transmitted and reflected 
modes on each side of the barrier at time iss (cf- Supple- 
mentary Note 4). In a realistic double- well potential, the 
wave packets undergo a complicated dynamics involving 
multiple transmissions and reflections, resembling Bragg 
beams splitters (see Fig. [4^l— b). Simulations calibrated 
with the measured trap parameters (see Supplementary 
Note 5 and Supplementary Method 3) predict a maxi- 
mum contrast of 80%, limited by the imperfect mode- 
matching in a double well with a finite barrier thickness. 
Taking into account interactions, no simple analytical 
model could be found. However, for weak interactions as 
in our case, numerical simulations yield a sine-like depen- 
dence of n(</>) and up to 70% contrast. Experimentally, 
we could achieve a contrast of 42 %. Possible reasons for 
this reduction are effects beyond mean-field and the cou- 
pling to the other spatial directions, both of which are 
not captured by our simulations. 

The steps described above are combined to perform the 
full interferometric sequence depicted in Fig. [T]d. Adjust- 
ing the phase shifter tilt e/h to —349 Hz, we record the 
final population imbalance z = n/N t as a function of 
the phase accumulation time t^. For each value of 
we typically perform 18 independent measurements of z, 
and compute the ensemble average and standard devia- 
tion of zftcf,). As expected, we observe that (z) oscillates 
at the frequency e/h. Phase diffusion manifests itself as 
a damping of the oscillations associated with increasing 
fluctuations of z. Note that the error bars of Fig. [TJ; rep- 
resent ±1 standard error of the mean. For > 20 ms, 
(z) remains close to because of the partial randomiza- 
tion of the phase while single experimental runs can still 
achieve ~ 40 % imbalance. 

To check the consistency of this data with that ob- 
tained by time-of-flight recombination (Fig. [3]) , we derive 
the expression for (z(t^,)} at the output of the recombiner 
from Eqs. (!]), Q and Q (cf. Supplementary Note 6). 
We find 

(z(t^) = Ce - A ^o/2 e -(^) 2 /2 sin( ^ + eU /h). (5) 

To compare this to the data of Fig. [T] we take the values 
of e, R and A0q obtained from the time-of-flight mea- 
surements (Fig. [3]). The contrast C = 0.42, on the other 
hand, is extracted from the data of Fig. [IJ by binning 
all the single-shot imbalances and fitting the histogram 
with the expected probability distribution function. We 
fit Eq. ([5]) to the data with ip and a constant offset as 
only free parameters (solid red line on Fig. [TJ;). The 
good agreement of the model with the measurement con- 
firms that the interrogation time of our interferometer is 
limited exclusively by phase diffusion. The black line in 
Fig. [Tfc is the signal expected in the absence of number 
squeezing. It underlines that the interrogation time is 
extended by the use of a number-squeezed state. 
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DISCUSSION 

We can now compare the phase estimation perfor- 
mance of our non-adiabatic recombiner to that of the 
usual time-of-fiight recombination scheme. An analysis 
of the measurements of Fig. [TJ: gives an intrinsic read-out 
noise of the non-adiabatic recombiner 8<p = 0.18 rad (cf. 
Supplementary Note 6), which is about twice the detec- 
tion noise of time-of-fiight recombination Scf)^ = 0.07 rad 
(see Methods). This difference is essentially due to the 
moderate contrast (42 %) of the non-adiabatic recom- 
biner. Were it increased to 100%, for instance through 
optimal control of the double-well potential [24] I53T 
155] ,the sensitivity of the non-adiabatic recombiner would 
compete with that of the time-of-flight technique. 

We also underline that phase estimation based on atom 
counting is more robust than phase extraction from inter- 
ference patterns. Firstly it does not require high imag- 
ing resolution, secondly, using a recombiner to convert 
the relative phase into a population imbalance allows the 
use of the precise atom counting techniques already avail- 
able [3T1I35JI371EIBEI]- Contrary to fringe fitting, it does 
not rely on time of flight, which makes it more suitable 
for a fully integrated device on a chip [55] ISH] • 

A more fundamental distinction between both methods 
is their potential sensitivity limits. Even though phase 
estimation based on a fit to the time-of-flight interfer- 
ence pattern can reach sub-shot-noise sensitivity, that is 
<50fit < <50sql oc 1/vAT, it is fundamentally bounded by 
gg This lower bound holds for any entangled 
state used in the interferometer. This non-trivial result 
underlines that the Heisenberg scaling of the phase sen- 
sitivity 5<fr oc 1/N will not be accessible to atom inter- 
ferometers using this phase estimation strategy. On the 
contrary, our non-adiabatic recombiner essentially works 
as a conventional optical beam splitter, and we thus con- 
jecture that it can reach the Heisenberg scaling [2"2l |2"4"] . 

In conclusion, we have developed a new intcrfcromc- 
try technique for trapped and spatially separated Bose- 
Einstein condensates based on a novel, non-adiabatic re- 
combiner. The geometry is reminiscent of the familiar op- 
tical Mach-Zehnder interferometer. Because its two arms 
are spatially separated, our device is sensitive to differen- 
tially applied forces. The tight confinement and precise 
positioning on an atom chip allows for local probing. We 
have demonstrated that using a number-squeezed state 
helps increasing the coherence time, which in turn is still 
limited to ~ 30 ms by dephasing due to interactions. 

Our observations illustrate a fundamental limitation of 
interfcrometry with trapped Bose-Einstein condensates 
for precision measurements: interactions dominate the 
randomization of the phase. This is in contrast to pre- 
vious experiments in which technical fluctuations of the 
energy difference or parasitic excitations during splitting 
were the dominant cause of dephasing [T7] [IS] . 



Strategies for creating optimal input states, for ex- 
ample by means of optimal control of the double 
well, will improve upon our present interferometer only 
marginally [23J . To go far beyond requires, in addition, 
phase estimation protocols based on Bayesian analysis 
of the full atom number distribution after the recom- 
biner [24. Developing these techniques could yield both 
a sub-shot-noise sensitivity and a long interrogation time, 
even in the presence of interactions. 

The non-adiabatic recombiner can also be seen as a 
new tool to perform rotations of the pseudo-spin describ- 
ing the many-body wave function [55j around the x axis 
of the Bloch sphere. Such a tool was missing for conden- 
sates trapped in double-well potentials. This may prove 
useful for the tomography of the many-body wave func- 
tion, e.g. to investigate strongly entangled states such as 
squeezed or Schrodinger cat states. 

METHODS 
Uncertainties and error bars 

The uncertainties given in the text, as in A0 = 0.16 ± 
0.01 rad or £^ = —7.8 ± 0.8 dB are statistical uncertain- 
ties calculated as follows: for directly measured quanti- 
ties they are the standard error of the mean assuming 
Gaussian distributions, for fitted parameters we use the 
standard error obtained from the fit, and for calculated 
values we use usual error propagation methods. 

Radio-frequency dressing 

The rf signal is sent through two wires which are par- 
allel to the main trapping wire [35] (cf. FigJlJ,). In order 
to turn the trap into a double-well potential, the fre- 
quency of the rf field is kept constant 30 kHz below the 
atomic Larmor frequency of 910 kHz due to the static 
fields (mj? = — > mp = ±1 transitions), while the am- 
plitudes in the two wires are ramped up to 52 mA peak 
to peak, creating an oscillating magnetic field of 1.7 G at 
the position of the atoms [TH] EI] ■ The distance between 
the wells can be tuned by changing the total rf intensity. 
To tilt the double well, the intensities in the rf wires are 
varied independently from each other. 

Time-of-fiight number and phase measurement 

At each point in the sequence, we can measure the 
relative phase and the atom number difference between 
the two condensates. To measure the phase, we abruptly 
switch off the trap, let the clouds overlap and image the 
interference pattern from below after 46 ms of free ex- 
pansion using our light-sheet fluorescence detection sys- 
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tem [37]. We extract the phase from the Fourier trans- 
form of the transverse density profile after integration 
along the longitudinal axis. Since the atoms are pro- 
jected onto the three Zeeman sub-levels of the F = 1 
manifold at switch-off, and to avoid the smearing of the 
interference patterns inevitably caused by stray fields, a 
magnetic gradient is applied during time of flight to sep- 
arate the clouds. The Fourier analysis is performed only 
on the atoms in rap = (more than half of the total 
atom number). 

Two different methods were used to access the atom 
number difference n = N-^ — ATr. For the data of Fig. [2^,, 
the dressing intensity was suddenly raised to increase the 
splitting distance and give opposite momenta to both 
clouds. For the data of Fig. ^Rp, the rf dressing intensity 
was quickly ramped down to initiate a transverse inward 
motion in the static harmonic potential. After the atoms 
acquired a sufficient velocity, the trap was switched off 
and the two clouds crossed and spatially separated. In 
both cases, we ended up with two well-separated clouds 
after time of flight, from which we could count the fluo- 
rescence photons in regions of equal size. 

Number distribution analysis 

To measure the atom number difference distribution, 
we performed the same analysis as in Ref . |44) . Assuming 
that the number of fluorescence photons detected from 
each atom is a random variable of mean p and variance 
a 2 , and including a random background of variance a 2 , 
the variance of the photon number difference between the 
two counting areas is 

As 2 =p 2 &N + Na 2 p + al (6) 

where £n = An/y/N is the number-squeezing factor. For 
our detector, the mean of the photon number distribution 
is p = 16 ± 1 photons/atom. Photon shot noise and the 
amplification noise of the EMCCD camera yield a 2 = 2p. 
Using more than 200 fluorescence pictures to measure 
As 2 and a\, we get the value £jy = 0.41 ± 0.04 (without 
correcting for imaging noise, we measure ^Ar jU nc = 0.52). 
Setting £n = in Eq. ^ enables to estimate the detec- 
tion noise (dotted red curve in Fig. [2^,). 

Phase distribution analysis 

To compare the phase distribution measured on Fig.[2]3 
with the one expected from a coherent state, we simu- 
lated the full phase extraction method. We first com- 
puted the one-body density distribution corresponding 
to the ballistic expansion of a symmetric superposition 
with 4> = 0, and generated the positions of the atoms for 
each realization by drawing AT to f independent atom po- 
sitions from this distribution. Since each atom is drawn 



independently, this procedure mimics a coherent state. 
We used AT to f = 500 to take into account the reduction of 
the number of atoms we extracted the phase from due to 
the Stern-Gerlach separation during time of flight. Tak- 
ing into account the photon shot noise, the amplification 
noise of the CCD camera and the optical resolution of our 
imaging system (we modelled the optical point-spread 
function by a Gaussian with RMS width of 10 /im, to be 
compared to the fringe spacing of 107 /im), we produced a 
large number (~ 1000) of interference patterns, on which 
we ran our Fourier phase extraction routine. We obtained 
a Gaussian phase distribution with standard deviation 
<5</>tof = 0.08 rad (dot-dash green curve of Fig. ^p), signif- 
icantly larger than the shot noise limit expected for an 
ideal detection system S^gQh = l/vA^ ~ 0.03 (dashed 
blue curve). 

A rough estimate of the phase detection noise is given 

1 /2 

by S<fia — (6<fi 2 of — <5^sql) — 0.07 rad, showing that 
detection noise dominates over projection noise. 



Test of uniform phase 

In Fig. [3)3 we used a Rayleigh non-uniformity test [45] 
to distinguish the points having a uniform phase distri- 
bution (black points) from the others (blue points). The 
blue points correspond to distributions whose p-values 
were below 0.05, indicating a probability less than 5% to 
be compatible with a uniform distribution. Only those 
points were fitted with Eq. ([2|. 
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SUPPLEMENTARY INFORMATION 

SUPPLEMENTARY NOTES 
Supplementary Note 1: adiabatic tilting 

The transfer from the horizontal trap to the tilted dou- 
ble well and the inverse operation are performed within 
3 ms each. Given the inter- well distance of 2 fim and the 
angles of ~ 10°, this corresponds to velocities at least 
30 times smaller than the widths of the velocity distribu- 
tions of the individual BECs, assuring that the transfer is 
sufficiently adiabatic to prevent radial dipole oscillations 
of the condensates. 



Supplementary Note 2: decoupling and phase 
diffusion 

From the trap calibration (see Supplementary Method 
1) we estimate a tunnel coupling energy |52| in the sym- 
metric double well just after splitting of the order of 
h x 0.1 Hz, high enough to prevent phase diffusion. How- 
ever, when we tilt the double well even by a few degrees, 
phase diffusion is observed. Since the currents in the 
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dressing wires are varied so as to keep the distance be- 
tween the minima constant, the coupling energy is not 
expected to vary much when the double well is tilted. 
Nevertheless, we checked that the phase diffusion rate 
we measured in the tilted trap was compatible with the 
one observed when splitting a condensate even more into 
a symmetric double well with vanishing tunnel coupling 
energy (< h x 10~ 8 Hz). This is consistent with simu- 
lations of the phase dynamics confirming that even for 
the lowest tilt angles we used, the energy detuning be- 
tween the wells was strong enough with respect to the 
tunnel coupling to effectively decouple the condensates 
and enable phase diffusion. 



Supplementary Note 3: Estimation of the chemical 
potential and phase diffusion rate after splitting 

Because the splitting is fast compared to the period Tb 
of the axial breathing mode (5 ms <C Tb = (v^^z) -1 — 
44 ms) , the process is not adiabatic for the axial motion 
but, on the contrary, can be considered sudden. Right af- 
ter splitting, the Thomas-Fermi (TF) radius is lager than 
its new equilibrium value because each BEC sees its den- 
sity suddenly divided by 2. Although this effect is par- 
tially compensated by a slight increase of v z and a change 
of the ID interaction constant during splitting, the ax- 
ial breathing mode is excited and the condensates are 
not at equilibrium. Consequently, their chemical poten- 
tials (e.g. entering Eq. ^) must be replaced by the "ef- 
fective" chemical potentials fi° s = Hdfa/dt, where <pi(t) 
is the time-dependent global phase of each BEC's wave 
function subject to a breathing mode (i € {R, L}, t = 
here corresponds to the end of splitting). With our pa- 
rameters, the BECs are well within the ID TF regime. 
Transversally, they are essentially in the ground state of 
the single-particle Hamiltonian (negligible interactions). 
The linear density profiles of the breathing BECs are well 
described by ID TF profiles whose radii Ri{i) = bi{t)R cq 
evolve according to the equation [61] 



h + ul bt = u 2 Jb 2 



(S7) 



Here, R° q denotes the equilibrium TF radii corresponding 
to each BEC. The corresponding phases do not evolve 
linearly with time but read 



&(t) = ^n(t)/h + eit/n, 



(S8) 



where Tj = f Q b~[ l (u) du, the /i° q are the chemical po- 
tentials that the condensates would have if they were at 
equilibrium, and the e, denote the potential energies at 
the center of the two wells. For short times, the phases 
can be approximated to first order in t/T^. Assuming 
bi(0) = 0, as expected for a sudden splitting we obtain 
(j)i(t) ~ i4 q t/(hbi(t = 0)), which yields 



ni(0) being the peak linear densities after splitting, and 
g the ID interaction constant after splitting. This ex- 
pression has a simple interpretation: one recognizes the 
usual expression of the chemical potential of an equi- 
librium quasi-BEC in the ID TF regime in which the 
equilibrium peak density has been replaced by the peak 
density right after splitting (one usually sets ej = 0). 

Symmetric splitting. Note that in the case of a sym- 
metric splitting we have nL.(0) = n^(0) = n(0)/2, where 
n(0) is the peak density before splitting and the two 
clouds breath in phase and with the same amplitude. 
The relative phase is simply given by Eq. (JlJ in which 
e = e R - e L . 

Imbalanced splitting. We now turn to the estimation of 



the phase diffusion rate. We can write Eq. (S9) as fi 



,cff 



gNin(0)/N because sudden splitting leads to identical 
density profiles before and after splitting, even if there is 
an imbalance (Ni ^ N/2). We readily obtain 



dfi 



M—N/1 



9- 



n(0) 



(S10) 



which we use to evaluate the phase diffusion rate in 
Eq. Q. From the knowledge of the trap frequencies, 
number of atoms, squeezing factor, and using a ID TF 
theory we get R — 46 ± 4mrad/ms. 

Another quantity of interest is the difference of chem- 
ical potential A/i = fi^ — /i^ which arises from a popu- 
lation imbalance z^O. From Eq. (S9) we obtain 



Ayu = gn(0)z = Nz 



N=N/2 



(Sll) 



.off 



gni(0) 



(S9) 



This enabled us to directly measure d/i/dAT. For this, 
the condensate was split in a tilted trap (e ^ 0) and then 
levelled back to the horizontal (e = 0). Adjusting the 
tilt angle for splitting gave control over the population 
imbalance which then drove the phase evolution in 
the levelled double well (see Fig. S5). We measured 
this effect for 6 values of z and indeed observed a 
linear dependence of A/i with z and found the slope 
gn(0)/h = 763 ± 53 Hz. This gives a third value of 
R = 57 ± 5mrad/ms, in agreement with both the 
expected theoretical value mentioned above and the one 
directly measured (see main text). 



Supplementary Note 4: Action of the recombiner 

Here we derive Eq. Q of the main text. In the trans- 
verse dimensions of our quasi-condensates |62j interac- 
tions can safely be neglected for the transverse spatial dy- 
namics, the non-linear term of the Gross-Pitaevskii equa- 
tion being negligible compared the the kinetic term. We 
assume that the trapping potential stays symmetric with 
respect to a vertical symmetry plane, V(—x,y,z,t) — 
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Supplementary Figure S5. Direct measurement of Afi(z) 
when the condensate is split to yield a population 
imbalance. The numbers labelling the points are the angles 
(in °) used during splitting to obtain the measured imbal- 
ance. The obtained population imbalance is responsible for 
Afj, measured in the levelled double well. The dashed line is 
a linea r fit t o the data from which gn(0) is extracted thanks 
to Eq. (Sll I. The horizontal error bars are ±1 standard error 
of the mean. The vertical error bars representing ±1 stan- 
dard error of the fitted energy difference are smaller than the 
points. 



V(x,y,z,t), which is expected from the atom chip sym- 
metry. Before recombination, the wave function de- 
scribing the two condensates has the same symmetry, 
except for a relative phase: iPbec(x, y, z, <bs = 0) = 
(ipL + iPr) A/2, with e^ipni-x^,?) = tph(x,y,z) = 
if)(x, y, z, £bs — 0)- From the symmetry of the poten- 
tial and propagating with the Schrodinger equation, one 
finds that the BEC wave function can be written at all 
time as 

V-bec = tyfo V, z, t) + e^i-x, y, z, t)] /y/2. (S12) 



The imbalance z = J x<0 \tpBEc\ 2 dr -_J^ >0 \ip B vc\ 2 dr 



can readily be calculated from Eq. (S12), which gives 
Eq. Q of the main text. The contrast is essentially 
given by the overlap of ip(x, y, z, t) and ip(—x, y, z, t). 
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Supplementary Figure S6. Distribution of population im- 
balances at all times. It exhibits the characteristic double- 
structure expected from sampling the sine of a uniformly dis- 
tributed phase. The contrast and the noise of the recom- 
biner are estimated from convolving the probability distribu- 
tion function of the sine of a random phase with a Gaussian 
noise. The error bars represent ±yki, where hi is the number 
of event in the bin of index i. 



patterns arising from the multiple reflections and trans- 
missions in the double-well potential. 

To validate these simulations, wc used the phase shifter 
to prepare states with 8 different relative phases and 
abruptly reduced the splitting similarly to the non- 
adiabatic recombination. After holding the atoms for an 
adjustable time tes, we eventually switched off the po- 
tential instead of separating the clouds. Since for our pa- 
rameters, the transverse density distribution after time 
of flight is almost homothetic to the initial transverse 
momentum distribution, this gives us access to the mo- 
mentum distribution at any time in the recombiner. We 
could compare it to the results of the simulation and 



found qualitative agreement for all phases (see Fig. S7 1 



Supplementary Note 6: Model for interferometric 
signal 



Supplementary Note 5: Simulations of the 
non-adiabatic recombiner 

We performed numerical simulations of the dynamics 
of the condensates during the non-adiabatic recombina- 
tion using a ID Gross-Pitaevskii solver in the double- 
well potential. The effective ID interaction constant was 
computed assuming a quasi-condensate profile longitu- 
dinally |63j and a Gaussian profile corresponding to the 
non-interacting ground state in the vertical transverse di- 
rection (see Fig. Ek-b). The simulations show complex 



To describe the output signal of the Mach-Zchnder in- 
terferometer (see Fig. [lp), we assume that the relative 
phase follows a normal distribution with time-dependent 
mean and standard deviation given respectively by 

\ A^) = y/Acfi + tRtt)*. ( * L6) 

The first equation describes the phase accumulation 
while the second accounts for phase diffusion. Assum- 
ing that the recombiner acts on the phase as 

n = N t C sin <f>, (S14) 
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binning all the single-shot population imbalances (see 
Fig. S6 ) and assuming that samples [0, 2n] uniformly. 



The probability density function of z is given by 



pdf (z) = 



1 



1 



(S16) 



for \z\ < C, and elsewhere. 

To account for technical noise of the Mach-Zchndcr 



interferometer, Eq. ( S16 1 is convolved with a Gaussian 
function of mean and standard deviation a z 64J. A fit 
to the data gives the value of C = 0.42 and a z = 0.075. 



Equation ( S15[ ) is used as a fit model for the mean 
population imbalance of Fig. [I]: with only free parameters 
ipo and a constant (slight) offset to account for imperfect 
balancing of the double well. Using the technical noise 
and the contrast extracted from the data of Fig. |S6| we 
estimate the noise on the phase estimation of the non- 
adiabatic recombiner as 5<j) = &n/\dn/ d<j)\$ = {) ~ o z jC ~ 
0.18 rad. 



SUPPLEMENTARY DISCUSSION 
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Supplementary Figure S7. Dynamics of the condensates 
in the recombiner. Using the phase-shifter, a state with 
a phase of 2.13 rad is prepared. Its dynamics in the recom- 
biner is monitored by time-of-flight imaging. From top to 
bottom: simulation of the transverse in situ density profile; 
simulation of the momentum distribution scaled to the 46ms 
time of flight, including the finite imaging resolution; mea- 
sured density profiles after expansion. The white dashed line 
corresponds to tss = 2.25 ms giving maximum contrast. We 
checked the agreement between data and simulation for 8 ini- 
tial values of the phase. 

we can integrate over the phase distribution to derive the 
expression for the expectation value of the normalized 
population imbalance z — (N-^ — N-r) /N t at any time: 

(z(t)) = Ce -A^/2 e -(fli,) 2 /2 sin ^ Q + et ^/ h y (S15 ) 

The initial phase spread is responsible for a reduced con- 
trast of the fringes at t = 0, while phase diffusion causes 
a Gaussian damping of the oscillations of (z). We set 
the phase diffusion rate R, the initial phase spread A</>§ 
and the detuning e to the values obtained from the in- 
terference pattern analysis (see Fig. [3]). The contrast is 
self-consistently extracted from the data of Fig. []J by 



Supplementary Discussion 1: ID effects 

It may seem surprising that phase fluctuations typical 
for quasi-condensates [?D] do not play a significant role 
here. Nevertheless we stress that the relevant quantity 
for interferometry is the relative phase between the con- 
densates and not the spatially fluctuating phase of a sin- 
gle BEC. Given our atom number, squeezing factor, and 
temperature, the relative phase correlation length [65] 
exceeds the length of the condensates. Experimentally, 
this is confirmed by the observation of a high contrast of 
the interference patterns integrated over the whole axial 
length of the cloud more than 80 ms after splitting (above 
60%). It therefore makes sense to define a single global 
relative phase between the two quasi-condensates and to 
neglect ID phase fluctuations. For longer times, however, 
the contrast of individual interference patterns starts to 
decrease, suggesting that ID effects might be an issue 
when trying to enhance further the interrogation time. 
This could be readily solved by using true condensates 
instead of quasi-condensates. 

Supplementary Discussion 2: Physical origin of the 
energy difference e 

The dominant contribution to the energy e is the differ- 
ence of gravitational potential between the two trap min- 
ima. A difference of 100 Hz corresponds to a vertical shift 
of 47 nm, which is just a fraction of the BEC's transverse 
width. Nevertheless, two other effects may contribute 
to e: different values of the magnetic field between the 
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two minima, investigated in Ref. [3T], and the variation 
of the trapping frequencies during the tilt, which shifts 
both the zero-point energies and the chemical potentials 
of the BECs. From simulations of the rf-dressed trap we 
estimate these two contributions to be below 10% and 2% 
respectively. Uncontrolled inhomogeneous electric fields 
close to the chip could also contribute to e. 

SUPPLEMENTARY METHODS 

Supplementary Method 1: Preparation of the 
condensate 

To minimize breathing oscillations during splitting, 
we prepare the initial quasi-condensate in a rf-dressed 
trap. The rf amplitude is slightly below the split- 
ting threshold i IT]. The potential is harmonic along 
the longitudinal (z) and vertical (y) directions, with 
trapping frequencies v z = 12.4 Hz and v y = 1.75 kHz 
respectively. In the horizontal, slightly anharmonic 
direction (x), simulations predict a spacing between the 
ground and first excited state of 1.02 kHz. By means 
of rf evaporative cooling, we prepare samples with less 
than 10% relative number fluctuations from shot to 
shot. The temperature is estimated from the density 
profile after time of flight using a stochastic model for 
phase-fluctuating quasi-condensates [66] . 



Supplementary Method 2: Splitting 

The double-well potential obtained after splitting is 
characterized by the trap frequencies of each individ- 
ual well v z = 13.2 Hz (longitudinal direction) and 
v y = 1.84 kHz (vertical transverse direction). Along 



the splitting direction, we measure for each well 
v x = 1.44 kHz. The two minima are separated by 2/im 
and simulations predict a barrier height of 3.7 kHz and a 
tunnel coupling energy of h x 0.1 Hz 52J. 

Supplementary Method 3: Calibration of the trap 
simulations 

We compute the geometry of the rf-dressed potential 
beyond the rotating-wave approximation by means of a 
Floquet analysis 167, 68 . The trap is first simulated from 
the knowledge of the chip design (see Fig. la of the main 
article) and the values of the control parameters. The 
latter are then adjusted in the simulation through inde- 
pendent measurements of the trap properties as follows: 
for the non-dressed trap we measure precisely the current 
in the trap wire and the current in the wires responsible 
for the longitudinal confinement. This leaves two free pa- 
rameters to fully define the static trap: the magnetic field 
Bq at the centre of the trap and the bias field B x , which 
are adjusted to match the measured atomic Larmor fre- 
quency and the measured transverse trap frequency. For 
the rf-dressed trap the absolute amplitudes of the rf cur- 
rent in each dressing wire are calibrated independently 
by measuring the trap bottom with one wire at a time. 
Then both wires are switched on and the relative phase 
between the rf currents is fine-tuned for the trap bottom 
of the dressed trap to match the simulation for different 
values of the global rf intensity. The good agreement of 
the simulation with the experiment is checked by mea- 
suring the distance between the wells (from the fringe 
spacing of the interference patterns) and the trap fre- 
quencies. Furthermore, measurements of the frequency 
of Josephson oscillations in strongly coupled double wells 
confirm the good calibration of the trap simulation used 
for example to simulate the recombiner. 



